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Low-lying Odd-parity States of the Nucleon in Lattice QCD 



o 

a; 



J3 
a; 

D 
43 



> 
O 

(N 

q 

ON' 

O! 
(N 



X 



M. Selim Mahbub, 1 Wasecm Kamlch, 1 Derek B. Leinwebcr, 1 Peter J. Moran, 1 and Anthony G. Williams 1 

(CSSM Lattice Collaboration) 

1 Special Research Centre for the Subatomic Structure of Matter, 
School of Chemistry & Physics, University of Adelaide, SA 5005, Australia 

The world's first examination of the odd-parity nucleon spectrum at light quark masses in 2 + 1 
flavor lattice QCD is presented. Configurations generated by the PACS-CS collaboration and made 
available through the ILDG are used, with the lightest pion mass at f 56 MeV. A novel method for 
tracking the individual energy eigenstates as the quark mass changes is introduced. The success 
of this approach reveals the flow of the states towards the physical masses. Using the correlation 
matrix method, the two lowest-energy states revealed are found to be in accord with the physical 
spectrum of Nature. 
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Lattice QCD is the only currently known ab-initio or 
first-principles approach to the fundamental quantum 
field theory governing the properties of hadrons, Quan- 
tum Chromodynamics (QCD). While the ground-state 
hadron spectrum of QCD is well understood, a deter- 
mination of the excited state energy spectrum is in the 
process of being revealed. Ultimately, the results can be 
compared with the existing experimental data and pro- 
vide predictions and motivations for future experiments. 

Hadron spectroscopy is dependent on the rich dynam- 
ics of the strong interaction. For example, the experi- 
mentally observed mass of the first positive-parity excita- 
tion of the nucleon, known as the Roper resonance, 
(1440) Pn, is surprising low compared to the lowest- 
lying negative-parity partner, (1535) Sn. This 
phenomenon is not observed in constituent or valence 
quark models where the lowest-lying odd-parity state oc- 
curs naturally below the first J p = i + excitation. 

Drawing on experimental results, we note that the 
Breit-Wigncr width of the iV _ (1535) state is « 150 
MeV, approximately half the width of the Roper 
Ar+(1440) Q. Furthermore, the branching fraction 
r(7riV)/r for (1535) j s 2/3 of the Roper. Together, 
these factors indicate a suppression of 1/3 in the cou- 
pling of 7rA^ to the AT - (1535) state relative to the Roper. 
Noting that the light ttN dressing makes the most im- 
portant self-energy contribution, it is anticipated that 
the self-energy dressings of irN for the AT - (1535) will be 
reduced to approximately 10% of that for the Roper. A 
consequence of this is to suppress the finite- volume ef- 
fects of the lattice QCD simulation which can otherwise 
lead to large energy shifts associated with the avoidance 
of energy-level crossings of the single and multi-particle 
scattering states. Similar arguments for the Af _ (1650) 
suggest irN self-energy contributions are suppressed to 
the 25% level. Thus, it is interesting to directly compare 
the results of our lattice QCD simulations with exper- 
iment and gain insight on the quark mass dependence 
of these states. While finite-volume effects are of resid- 



ual interest in this investigation, understanding the finite 
volume effects on these states and linking them to the res- 
onances of Nature is a long term program of the lattice 
QCD community. 

The experimentally observed nearly-degenerate Sn 
(1535) and (1650) states are in agreement with the sim- 
ple quark-model predictions based on SU(6) symmetry. 
Therefore, looking at the low-lying energy states 

and their structure from the first principles approach is 
potentially very revealing. Some recent full QCD results 
can be seen in Refs. [2h7[. Herein, it will be interest- 
ing to explore the physics associated with the dynamical 
fcrmion loops of full QCD, this time at very light quark 
masses. 



The correlation functions for the 



states are short- 

lived compared to the lighter N± ground state. Thus 
it is important to adopt a method that can isolate the 
effective-mass plateaus at early Euclidean times. The 
variational method [1, Q is the state-of-the-art approach 
for achieving this in lattice hadron-spectroscopy calcula- 
tions and is adopted here. Through a generalized eigen- 
value analysis of a matrix of correlation functions, the 
process enables one to create highly optimized interpo- 
lating fields designed to excite a single energy eigenstate 
of the QCD Hamiltonian. The masses of the energy states 
are then obtained through a standard effective-mass anal- 
ysis flfl ] providing a robust approach for extracting the 
energy states at early Euclidean times. 

In this paper, we utilize the established approach of 
Refs. d, 0] to explore the low-lying energy states 

in full QCD. In doing so, a novel method has been devel- 
oped to track the energy eigenstates from heavy to light 
quark masses. The method is particularly useful when 
the energy-states are nearly degenerate. 

The two-point correlation-function matrix for p = 



2 



can be written as 

G « (*) = E Tr ^p i r ± < 1 x,(o) 1 )}, (1) 

= ^AfA« e -™^ (2) 

a. 

where Dirac indices are implicit, Xf and A" are the cou- 
plings of interpolators Xi an d Xj a t the sink and source 
respectively, a enumerates the energy eigenstates with 
mass m a , and r± = (70 ± l)/2 projects the parity of 
the eigenstates. A linear superposition of interpolators 
4> a = Xj u< j creating state a provides the relationship 



Gijih + At) u° = er m ^ G I3 (t ) u? , 



(3) 



from which right and left eigenvalue equations are ob- 
tained 



[(G(to)r 1 G(t + At)]« u? = c° 
v? [G(t + At) (G(t„))- 



= c<*vf, 



(4) 
(5) 



with c a = e~ m " . The vectors u" and vf diagonalizc 
the correlation matrix at time to and to + At making 
the projected correlation matrix, v"G^(t)Uj oc S a/3 . The 
parity and eigenstate projected correlator, 



(6) 



is then analyzed to obtain masses of energy-states. 

A eigenvector analysis of a symmetric matrix having 
orthogonal eigenvectors can be constructed by insert- 
ing G(t o y* G(t ) 3 = I, in Eq. (gj) and multiplying by 
G(to)? from the left, 



G(t ) *G(t Q + At)G(t r*G(t )*u a 
G(t Q )~% G(t + At) G(t ) 



c a G(t yu a 7 

(7) 



w = c w 



(8) 



where, w a = G(t ) 3 «" and [G(t )~ 2 G(t + 
— - 

At)G(to) 2 ] is a real symmetric matrix, with orthog- 
onal eigenvectors w a . The vector u a may be recovered 

_ 1 

from the w a via u a = G(to) 2 w a . 

The PACS-CS 2 + 1 flavor dynamical-fermion con- 
figurations [ll| made available through the ILDG [12 1 
are used herein. These configurations use the non- 
perturbatively C(a)-improved Wilson fermion action and 
the Iwasaki-gauge action [l3|. The lattice volume is 
32 3 x 64, with j3 = 1.90 providing a lattice spacing of 
a = 0.0907 fm and a physical volume of w (2.90 fm) 3 . 
Five values of the (degenerate) up and down quark 
masses are considered, with hopping parameter values 
of K ud = 0.13700,0.13727,0.13754,0.13770 and 0.13781, 
corresponding to pion masses of m n = 0.702, 0.572, 



S2 











1 








I 

I 




I 


I 


i 




I 








i 






i 


i 




z 


i 




i 




i 


I 




• 


■ 


■ 


m 


1 


■ 


■ 










■ 


. 









17 



18 
t„ 



19 



FIG. 1. (Color online). iVi energy-states from a 4 x 4 cor- 
relation matrix analysis of the xi interpolator at the lightest 
pion mass of = 156 MeV. The variational parameters to 
and At are shown at the major and minor tick marks respec- 
tively. The LHS of the diagram refers to to — 17, while the 
RHS is for t = 18. 



K s = 0.13640. Gauge-invariant Gaussian smearing [14J is 
used at the fermion source and sink with a fixed smearing 
fraction and four different smearing levels including 16, 
35, 100, and 200 sweeps 0,0] ■ 

The complete set of local interpolating fields with three 
different spin- flavor combinations for the spin-i nucleon 
are considered herein, 



Xi(x) =e abc (u Ta (x)C l5 d b (x))u c (x) 



X2(x) 

Xi{x) 



abc f n Ta 



(u la (x)Cd°(x)) l5 u c (x), 



(u la (x) C l5 -/ 4 d b (x))u c (x). 



(9) 
(10) 
(11) 



0.413, 0.293, 0.156 GeV [ll|; for the strange quark 



Each interpolator has a unique Dirac structure giving rise 
to different spin-flavor combinations. Moreover, as each 
spinor has upper and lower components, with the lower 
components containing an implicit derivative, different 
combinations of zero, one and two-derivative interpola- 
tors are provided. The interpolator \i is the time compo- 
nent of the local spin-| isospin-i interpolator which also 
couples to spin-i states. It provides a different linear 
combination of zero- and two-derivative terms comple- 
mentary to xi- 

In Fig. [TJ projected masses of the states are pre- 

sented from a 4 x 4 correlation matrix constructed from 
the interpolator xi an d four different smearing levels. 
The dependence of the results on the variational param- 
eters to and At is illustrated. While the lowest energy- 
state is almost independent of to and At, the excited 
states show some dependence at smaller to an d At val- 
ues. The energy-states at (to, At) = (18,2) provide the 
best balance between the systematic and statistical un- 
certainties [f| and these parameters are therefore selected 
for our numerical study. 

In Fig. [5] we show results for the lowest energy-state 
from dynamical and quenched (l5| QCD simulations. As 
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FIG. 2. (Color online). Dynamical and quenched results for 
the lowest Ni energy-state using the Xi interpolator. 
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FIG. 3. (Color online). N\ energy-states from an 8 x 8 cor- 
relation matrix of \\ and X2 interpolators, for the pion mass 
range of 156 to 702 MeV. The physical spectrum [l[ is 

shown at the far left. 



anticipated, the quenched and dynamical results are in 
agreement in the heavy quark-mass region. However, in 
the light quark-mass regime the results are significantly 
different as the effects of the light sea quarks become 
increasingly important. Only the dynamical results ap- 
proach the physical value and this provides strong evi- 
dence for the non-trivial role of light sea-quark degrees 
of freedom to the structure of Nature's hadron spectrum. 

To explore the nearby second energy state, Sn (1650), 
we extend our analysis to include the interpolators X2 an d 
X4 with a variety of smearing levels. The results of an 8 x 
8 correlation-matrix analysis of xi and xi interpolators 
with four levels of smearing arc presented in Fig. [3l 

The flow of the lowest two energy states towards the 
physical values is remarkable. The results at the two 
heaviest pion masses sit close to the scattering S-wave 
N+ir threshold indicating that the results for these heav- 
iest pion masses may be scattering states. However, in 
the light quark-mass region these states move above the 
threshold. 



It is likely that the three-quark interpolators used 
herein have relatively small couplings to the scattering 
states at the light quark masses relative to the states 
observed and are not resolved in the correlation-matrix 
analysis. Hence, a combination of five- and three-quark 
interpolators in a correlation matrix analysis is highly de- 
sirable for future investigations to ensure better overlap 
with the multi-particle states. This type of novel work 
using the stochastic LapH method is in progress [l6[ . It 
is because the coupling to multi-particle states at light 
quark masses is heavily suppressed, that it is meaning- 
ful to compare lattice results with the central values of 
experimentally measured hadron resonances. 

A similar situation prevails for the second pair of states 
in the spectrum, where the states sit close to the P-wave 
En + + M n and E„ + E n + Mm threshold scattering 
states with back-to-back momenta of one lattice unit, p = 
(2n/L x ,Q,0). The apparent flow of these states in the 
light-quark region toward the physical Sn (2090) state is 
also interesting. 

In presenting the results of Fig. [3] and assigning sym- 
bols to each of the energy levels observed at a particular 
quark mass, it is necessary to track the evolution of the 
states from one quark mass to the next. We have done 
this through a consideration of the evolution of the eigen- 
vectors as the quark mass is changed. 

Consider M interpolating fields making an M x M 
parity-projected correlation matrix G(t) and its as- 
sociated symmetric generalized eigenvalue equation of 
Eq. (|8|). Using the normalization \ w ?\ 2 = 1> tnc 
quantity w a (m q ) ■ w^(m q ) = S a p. This feature enables 
the use of the generalized measure 

W aP (m q ,m q ,) = w a {m q ) ■ w p (m q> ) 



(12) 



to identify the states most closely related as we move 
from quark mass m q to an adjacent quark mass m q > . The 
state numbers a and /? are assigned in order of increas- 
ing projected eigenstate energy at the quark masses m q 
and m q > respectively. Typical results for this generalized 
measure of eigenvector overlap are presented in Table |TJ 

For each value of state index a there is only one value of 
(3 where the magnitude of the entry is significantly larger 
than all others and approaching unity. The most relevant 
entries for consideration are the immediate neighbors of 
a where a crossing of the eigenvectors moves the largest 
entry off the diagonal. 

This measure provides a clear identification of how 
states in the spectrum at quark mass m q are associated 
with states at the next value of quark mass, m q i . For ex- 
ample, the results of Table |T] indicate the first four states 
at m q i appear with the same ordering in the spectrum as 
observed at m q , the fifth state at m q ' is associated with 
the sixth state at m q and vice versa and similarly for the 
seventh and eighth states. We note that while the central 
values of the energies have changed order, the error bars 
are sufficiently large that one cannot conclude that an 
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TABLE I. The scalar product w a (m q ) ■ w 8 (m q i) for «: = 
0.13754 (m„ = 413 MeV) and k' = 0.13770 (m, = 293 MeV) 
for an 8 x 8 correlation matrix of xi and X2 with four different 
levels of smearing. State numbers a and /3 correspond to row 
and column number, respectively. 
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0.40 
-0.01 
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0.04 
0.01 
0.00 
0.01 
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-0.01 
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FIG. 4. (Color online). The components of the eigenvector 
u a providing the amplitude for each interpolating field at the 
source for creating the state a. The states are labeled by the 
eigenvector (EVect) number with the ordering as provided 
in Fig. [3] at the heaviest quark mass. For each EVect, the 
eigenvector components are plotted from left to right in order 
of increasing quark mass. In the legend, (ui, U2), («3, U4), 
(us, uq) and (117, us) correspond to the smearing-sweep lev- 
els of 16, 35, 100 and 200 respectively. Odd numbers in the 
subscripts correspond to the contribution from the xi inter- 
polator, whereas, even numbers correspond to X2- 



avoided energy level crossing has taken place in moving 
from quark mass m q to m q > . 

The components of the eigenvector u a , providing the 
amplitude for each interpolating field at the source for 
creating the state a, are provided in Fig. [4j A non- 
trivial contribution from both the xi and xi interpolators 
for the lowest two energy-states is evident. The scalar- 
diquark interpolator xi dominates the lowest energy- 
state. On the other hand, both xi and xi interpolators 
have large contributions to the second energy state where 
their strengths appear with opposite signs. The eigenvec- 
tor components typically display a slow evolution as the 
quark mass is changed. 

The energy-states for our complete analysis are pre- 
sented in Fig. [5] The results are drawn from two 8x8 
correlation-matrix analyses for pairs of xi, Xi and xij Xi- 
The matrices are formed with each interpolator having 
four levels of smearing. Whereas the xi; Xi and X2, Xi 
analyses reveal a similar spectrum, four new states are 
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FIG. 5. (Color online). Masses of 12 low-lying JVi energy 
states from two 8x8 correlation matrices of Xi, X2 an d Xi> Xi- 
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FIG. 6. (Color online). The quark mass dependence of the 
lowest two lowest-lying states are compared with the 

S-wave scattering threshold. 



revealed in the xij Xi analysis providing the resolution 
of 12 low-lying states in our analysis. 

In the quark model based on SU(6) spin-flavor sym- 
metry, the odd parity (1535) and (1650) states belong 
to the negative parity, L = 1, 70-plet representation of 
SU(6). As three spin-i quarks may combine to a total 
spin of s = I or | , the L = 1 state can couple two differ- 
ent ways to provide a J = | state, hence providing two 
orthogonal spin-i states in the L = 1, 70-plet represen- 
tation. Both of these states have a width of « 150 MeV. 
The lowest two energy states revealed here are similarly 
close in mass, as illustrated in Fig. [5l in accord with the 
SU(6) quark model. 

These two lowest-lying states are presented in 

Fig. [5] in comparison with the 5-wave scattering thresh- 
old. These lattice results, providing the first examination 
of the odd-parity nucleon spectrum at a pion mass as low 
as 156 MeV, display remarkable agreement with the phys- 
ical values. They represent a significant achievement for 
lattice regularized QCD in describing Nature. 

Although both these low-lying states are quite similar 
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at the two heaviest quark masses, their approach to the 
/sica] values in the light qna.rk-ma.ss region are Hiffer- 



[4] G. P. Engel, C. B. Lang, M. Limmer, D. Mohler. and 
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ent . I Significant chira.l curvature is evident in particular 



A. Schafer (BGR [B ern-Graz-Regensburgl), |Phys. Rev 



for the second state. It will be interesting to explore the 
mass dependence of these states using effective field the- 
ory techniques and to repeat these studies on matched 
lattices of different volume when they become available. 
Future studies will endeavor to observe the multi-particle 
scattering states and determine the resonance parameters 
of these states from the first principles of QCD. 
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